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Abstract. Experimentally assessing the preferred body temperature (T „„;) of ectothermic animals is important to under- 
stand thermal adaptation. In lizards, this variable is usually estimated by measuring body temperature in thermal gra- 
dients. To quantify the extent to which different experimental setups influence the inferred T „values we submitted 65 
individual lizards of three species to randomized tests using six different experimental setups at Oukaimeden, Morocco, 
including setup variants similar to those that have been most often used in lacertids. Among-treatment differences were 
substantial. Using an infrared bulb as heat source in combination with artificial cold lighting yielded about 5°C lower Tç 
estimates than photothermal treatments with an incandescent bulb as heat and light source, possibly because lizards ther- 
moregulated differently without a visual cue related to the heat source, or due to the absence of a natural photoperiod. 
Photothermal assays in which T „was assessed by hourly cloacal measurements over a 10h period yielded 2.2°C lower 
T „estimates than 2h treatments where body temperature was measured every minute with a thermocouple attached to 
the belly. This probably reflects that the 2h treatments targeted lizards in the initial warming-up phase, whereas the 10h 
treatment attempts to capture the preferred temperature of a lizard over its entire daily cycle including phases of inactiv- 
ity. Lastly, we observed large differences among treatments with contact thermometers versus infrared laser thermometer 
measurements, calling for caution when the latter are used with artificial heat sources. Our data do not provide thorough 
tests of the physical, behavioural or physiological causes underlying the observed differences between treatments, but illus- 
trate that for meta-analyses where detailed comparisons are needed, a rigorous consideration of the optimal experimental 
setup and its consistent use will remain necessary. 


Key words. Squamata, Lacertidae, Podarcis, Atlantolacerta, Scelarcis, Morocco, thermoregulation, preferred temperature, 
field body temperature. 


Introduction ecology and conservation, and depends on robust experi- 

mental methods. This applies for instance to lizards, which 

Assessing the thermal preferences of ectothermic animals have served as a model group to develop and test the hy- 
is crucial to numerous questions in evolutionary biology, pothesis of faster adaptive processes to cold than to heat 
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(Aravjo et al. 2013, MuÑoz et al. 2014), because their ca- 
pacity for thermoregulatory behaviour may serve as a buf- 
fer in the case of cooling. Macroecological studies have 
elucidated that environmental temperatures predict lizard 
life history better than their body temperatures (MEIRI et 
al. 2013). However, in one recent study, estimations of oc- 
currence based on activity outperformed those based on 
temperature alone, which suggests other environmental 
factors (e.g. humidity) are also involved (CAETANO et al. 
2020). Among lizards, forest-dwelling species are mostly 
thermoconformers and cold-adapted species are usually 
thermoregulators, and both groups might suffer from ther- 
mal and/or hydric stress, and undergo declines, with glob- 
al warming and aridification (SINERVO et al. 2010, DIELE- 
VIEGAS et al. 2018, GARCIA-PorTA et al. 2019). Analyzing 
these ecological, evolutionary and conservational ques- 
tions depends on inferring a series of physiological traits 
that inform on the interactions of lizards with their ther- 
mal environment (CLUSELLA-TRULLAS & CHOWN 2014). 
These encompass thermal limits for survival (critical tem- 
peratures: CT a CT.) but especially for activity (volun- 
tary thermal maximum and minimum: VT a VThi CA- 
MACHO et al. 2018), field body temperatures measured in 
active lizards in the wild (T,), optimal performance tem- 
perature (i ) and preferred body temperature (Ty ap: Es- 
pecially this latter variable, T „p has played a prominent 
role in understanding thermal adaptations (reviewed in 
CLUSELLA-TRULLAS & CHOWN 2014), and is a crucial com- 
ponent of models of climate change responses of lizards 
(e.g., SINERVO et al. 2010, GARCÍA-MUÑOZ & CARRETERO 
2013). In many lizards, T „is also highly correlated with 
other physiological optima “such as running speed or diges- 
tion (BAUWENS et al. 1995). 

Experimentally, T,_,is measured by exposing lizards to a 
temperature gradient i in the supposed absence of costs and 
constraints of operative environmental temperatures that 
are present in field conditions (BAKKEN et al. 1985; but see 
SANNOLO & CARRETERO 2019 for hydric costs), and record- 
ing their body temperature over time (CAMACHO & RUSCH 
2017). Traditionally, lizard body temperature is measured 
as cloacal temperature (HUEY & WEBSTER 1975, 1976, HUEY 
& PIANKA 1977), an approach also used in T „experiments 
by many researchers (e.g., CARRETERO & LLORENTE 1995, 
MARQUEZ et al. 1997, VERÍSSIMO & CARRETERO 2009, CAR- 
RETERO 2012, OSOJNIK et al. 2013, ZAMORA-CAMACHO et al. 
2014, Kapsatas et al. 2016). Other studies monitored body 
temperatures with a thermocouple attached to the outside 
of the body (e.g., CALSBEEK & SINERVO 2002, AGUADO & 
BRAÑA 2014, MORENO AZOCAR et al. 2016, GARCIA-PORTA 
et al. 2019), implanted into the body cavity or underneath 
the skin (e.g., LICHT 1968, COWGELL & UNDERWOOD 1979), 
or measured temperature of the body surface using a laser 
thermometer (CARRETERO 2012, BOUAZZA et al. 2016). Re- 
markably, CARRETERO (2012) found that lizard body tem- 
peratures measured by infrared vs. contact thermometers 
were only poorly correlated and with a biased relation. In 
contrast, BARROSO et al. (2016) obtained a good correlation 
and minimum bias in the comparison of thermographic 
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images with cloacal temperature measured by contact 
thermometers. As a further potentially confounding fac- 
tor, even in small lizards regional heterothermy occurs, i.e., 
different body parts can have different temperatures (SAN- 
NOLO et al. 2014, BARROSO et al. 2016, 2020). 

The kind of thermal gradient used also differs substan- 
tially between studies (reviewed in CAMACHO & RUSCH 
2017). Some researchers use incandescent light bulbs as 
heat source to create a photothermal gradient (e.g., LICHT 
1968, GVOZDIK 2002, AGUADO & BRAÑA 2014, KIRCHHOF 
et al. 2017, GARCIA-PorTa et al. 2019), while others ad- 
vocate the use of uniformly distributed (natural) lighting 
originating from windows, and use infrared (IR) bulbs to 
create a purely thermal gradient (e.g., BOWKER and JoHN- 
SON 1980, ANGILLETTA 2001, LAILVAUX et al. 2003, VERÍS- 
SIMO & CARRETERO 2009, CARNEIRO et al. 2015). It is likely 
that such differences between heat and light sources lead to 
behavioural differences of the lizards in the gradient, mir- 
roring differences in thermoregulatory behaviour in shad- 
ed vs. exposed environments (e.g., HUEY 1974). 

A further difference between studies refers to the time 
during which lizards are exposed to a temperature gradi- 
ent. Individuals can be kept full days up to weeks in the 
gradient with occasional temperature measurements 
(which often involves catching the lizard out of the gra- 
dient thereby influencing its behaviour and body temper- 
ature via convection and removal from the original tem- 
perature source) (LICHT et al. 1966, LICHT 1968, SIEVERT 
& HUTCHISON 1988, 1989, CARNEIRO et al. 2015, ORTEGA 
et al. 2016). Alternatively, lizard body temperature is mea- 
sured over a relatively short period of a few hours, ran- 
domized throughout the day, mostly with continuous tem- 
perature measurements, to minimize the possible effect of 
diel variation in ine op (€-g-» PARANJPE et al. 2013, YANG et al. 
2008, Lī et al. 2009, WANG et al. 2013, KIRCHHOF et al. 2017, 
GarciA-Porta et al. 2019). 

Given the importance of cross-taxon comparisons and 
meta-analyses to understand evolutionary and macroeco- 
logical patterns of lizard thermal physiology (SINERVO et al. 
2010; MEIRI et al. 2013, DIELE- VIEGAS et al. 2018, GARCIA- 
Porta et al. 2019), it is of high importance to assess how 
these methodological differences influence the inferred 
values of T „p While the reviews of CLUSELLA-TRULLAS & 
CHOWN (2014) and CAMACHO & RuscH (2017) highlighted 
the wide variation of methods used, a rigorous experimen- 
tal approach is necessary to quantitatively assess impacts of 
using different methods (e.g., CARRETERO 2012, BARROSO 
et al. 2016). 

Thus far, methodological studies on lizard thermal phys- 
iology have (i) demonstrated differences in thermal pref- 
erence at different times of the day (e.g. COWGELL & UN- 
DERWOOD 1979, RISMILLER & HELDMAIER 1982, SIEVERT & 
HUTCHINSON 1988, FIRTH & BELAN 1998), (ii) shown differ- 
ences in preferred temperature depending on the relative 
position of heat vs. light source in diurnal and nocturnal 
lizards (SIEVERT & HUTCHINSON 1988, 1989), (iii) compared 
measurements of infrared laser vs. contact thermometers 
(CARRETERO 2012) and thermographic images vs. contact 
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thermometers (BARROSO et al. 2016), and (iv) compared 
temperature measurements at different lizard body parts 
(SANNOLO et al. 2014, BARROSO et al. 2016, 2020). However, 
experimental comparisons of approaches to T „inference 
from long vs. short segments of the diel activity period, 
or from trials with infrared vs. incandescent bulbs as heat 
source, have not been carried out yet. Ideally, such compar- 
isons should be carried out in the same study period and 
lizard species, and if possible the same lizard individuals. 

Here, we empirically quantify differences among T 
measurements obtained from different experimental ap- 
proaches on the same organisms. Our study focuses on 
three species of the family Lacertidae, an Old World clade 
that has been intensively explored for thermal ecology. 
CARRETERO (2012) and GARCIA-PORTA et al. (2019) have 
highlighted differences in T, and T,,, values of lacertids 
obtained using different experimental approaches. We ex- 
pand these data by assessing T „of the same individual liz- 
ards in six different settings reflecting the most common- 
ly used approaches in these lizards, mainly to understand 
the impact of (a) full day hourly measurements of cloacal 
temperature vs. short-term measurements by ventrally at- 
tached thermocouples, and (b) usage of a photothermal 
gradient vs. thermal gradient with uniform lighting; and 
furthermore, (c) the difference between contact and la- 
ser thermometers, as previously assessed by CARRETERO 
(2012). The main goal of our study is to describe the differ- 
ences between various experimental approaches that have 
often been applied to estimate T „in lacertids and/or other 
lizards. We provide hypotheses that could explain the ob- 
served differences among methods but do not aim at pro- 
viding functional explanations within a thorough hypoth- 
esis-testing framework. 


Methods 


Experiments were carried out in Oukaimeden, Morocco 
(approximate geographical coordinates 31.2075°, -7.854°, 
ca. 2600 m above sea level), with lizards of three sympat- 
ric species: Atlantolacerta andreanskyi (WERNER, 1929), 
Podarcis vaucheri (BOULENGER, 1905), and Scelarcis perspi- 
cillata (DUMERIL & BIBRON, 1839). Active lizards were di- 
rectly captured in the field in the breeding season, April 
and May, corresponding to the lizards’ main activity peri- 
od. Overall, 30 A. andreanskyi, 23 P. vaucheri, and 12 S. per- 
spicillata, all adult males to exclude biases due to sex, re- 
production or body condition in thermal preferences, were 
included in the experiments (S’KHIFA et al. 2020). Active 
specimens were captured in the field on sunny days, op- 
portunistically during random searches between 9-17 h lo- 
cal time. Field body temperatures were measured immedi- 
ately upon capture as (i) cloacal temperatures as described 
in treatments 3 and 4 below, and (ii) as dorsal surface tem- 
peratures as described in treatments 5 and 6 below. 

Our methodological approach relies on determining se- 
lected body temperature in a thermal gradient as a proxy 
for preferred body temperature (T). All experiments 


were conducted in an experimental setup consisting of five 
individual tracks, each of which ca. 120 cm in length and 20 
cm in width, corresponding to a surface of ca. 2400 cm’. 
Experiments started at 8:30 h in the morning, and again at 
14:30 h in the afternoon for the 2 h treatments, to mirror the 
lizards’ regular activity period. Before and after the experi- 
ments, lizards were kept for several hours at an air tempera- 
ture of 10-15°C (corresponding to natural air temperatures 
outside of the sun), providing water and food ad libitum. 
In all treatments, the gradient ranged from surface and air 
temperatures of ca. 10-15°C in its coldest part to very hot 
temperatures (> 40°C) in the immediacy of the (infrared 
or incandescent) bulbs used as heating source, and this also 
approximately represents the minimum and maximum 
potential operative temperatures that these small lizards 
could attain in the gradient if spending sufficient time in 
the extreme parts of the gradient. After completion of the 
experiments, lizards were released at the site of initial cap- 
ture. When lizards showed signs of fatigue or non-natural 
behaviour, or when completion of all treatments was lo- 
gistically not feasible within five days, they were released 
back at the capture sites. Each lizard was subjected to each 
of the six treatments following a randomized order. Alto- 
gether, 65 lizards were studied, and for 21 A. andreanskyi, 
20 P. vaucheri, and 11 S. perspicillata complete data sets for 
each of the following six treatments were obtained: 

(1) photothermal-zh: an incandescent 100 W light bulb 
(full spectrum) was suspended 30 cm above one end of a 
track, to create a photothermal gradient of approximate- 
ly 15 to 55°C at ground level. Body temperature was deter- 
mined every minute by ultra-thin T-type thermocouples 
(OMEGA 5SCTT-T-40-72, diameter = 0.076 mm, Nor- 
walk, Connecticut, USA) affixed with medical tape to the 
lizards’ venter and connected to an 8-Channel USB Ther- 
mocouple Data Acquisition Module (OMEGA TC-08; res- 
olution < 0.1°C). Lizards were allowed unrestricted move- 
ment within their individual gradient. Lizard body tem- 
perature was measured over a period of up to 2.5 hours. 
Lizards were randomly assigned to either morning or after- 
noon treatments. The initial 30 minutes of the experiment 
were discarded as acclimation time, resulting in a maxi- 
mum of 120 individual measurements (= 2 hours). Indi- 
viduals that remained inactive at the cold end of the gradi- 
ent for more than 50 minutes were not included in the final 
analysis assuming they were not thermoregulating. Tem- 
porarily entangled or detached thermocouples were fixed 
and reattached, the respective data points were discarded, 
and the experiment was continued. This treatment corre- 
sponds to the one used in GARcIA-Porta et al. (2019). 

(2) thermal-2h: as in the previous treatment, but instead 
of an incandescent light bulb, a 75 W IR bulb was used and, 
separate, uniform cold lighting provided for the entire gra- 
dient and for the entire time of the experiment. 

(3) photothermal-1oh: as in photothermal-2h (incan- 
descent bulb), but experiment was run for 10 hours, lizards 
caught once per hour, and their cloacal temperature mea- 
sured within 20 s with a digital quick-reading thermome- 
ter (GHM-Greisinger: GTH 1170, GHM, Regenstauf, Ger- 
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many), introducing the thin external probe (K-type thermo- 
couple) ca. 2 mm into the lizard’s cloaca within 10 s after cap- 
ture. Temperatures were recorded with a resolution of 0.1°C. 

(4) thermal-1oh: as in photothermal-1o h, but with a 
75 W IR bulb and uniform cold lighting. This treatment 
largely corresponds to the one extensively used in previous 
studies of lacertid thermophysiology (e.g., CARRETERO et 
al. 2006, VERÍSSIMO & CARRETERO 2009, GARCIA-MUNOZ 
& CARRETERO 2013) except for the use of artificial rather 
than natural light as in the quoted studies, and narrower 
gradient setups (30 cm in the quoted studies). 

(5) photothermal-1ohL: as in photothermal-ioh, but 
temperatures measured each hour with an infrared laser 
thermometer (Fluke 62; Fluke Europe B.V., Eindhoven, 
The Netherlands) dorsally near the tail base in line with 
the body axis from approximately 300 mm distance (Car- 
RETERO 2012, AGUADO & BRANA 2014, LARA-RESENDIZ et 
al. 2015). The Fluke 62 measures with a resolution of 0.1°C 
and an accuracy of +1.0°C or +1.0% of reading with a sen- 
sor diameter of 38 mm (at 300 mm distance). 

(6) thermal-1ohL: as in photothermal-1oh, but with an 
IR bulb and uniform lighting, and temperatures measured 
each hour with an infrared laser thermometer as in photo- 
thermal-1ohL. 

In the 2h treatments, the lizards were observed during 
the entire time of the experiment, and all time periods in 
which thermocouples were detached or entangled were 
noted, and later excluded from analysis. We also excluded 
all obvious outliers corresponding to erroneous measure- 
ments from the raw data set. We then calculated individual 
medians, yielding one datapoint per individual and treat- 
ment (i.e., typically six data points per individual). 

As usual in these kinds of data, they were skewed to- 
wards higher values and therefore not normally distrib- 
uted (Shapiro-Wilk tests, p > 0.05 in comparisons of the 
full data set and subsets by method), even after log-trans- 
formation. We therefore used Poisson-distributed Gener- 
alized Linear Mixed Models (GLMMs) which are robust 
against non-normally distributed data, as implemented in 
JMP 13.0 (SAS Institute) with the GLMM add-in (https:// 
community.jmp.com/ts/JMP-Add-Ins/Generalized-Lin- 
ear- Mixed-Model-Add-in/ta-p/284627) to assess the in- 
fluences of various predictors on inferred T „p Specifically, 
the following models were calculated: 

Model 1: To understand if time period (morning or af- 
ternoon) played a role in explaining T „inferred from the 
photothermal-2h and thermal-2h treatments, an initial 
GLMM was conducted with T ¿as response variable, ex- 
perimental method (treatment), s species and time period as 
predictor variables, and lizard individual as random factor 
(no interactions specified), only including data from the 
two 2h treatments. A total of N=106 Ta , Values were in- 
cluded in this model. 

Model 2: A GLMM with pi ¿as response variable, ex- 
perimental method and species as predictor variables, liz- 
ard individual as random factor (no interactions specified), 
for all treatments. A total of N = 346 Ti , Values were in- 
cluded in this model. 
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Model 3: A variant of Model 2 with the same factors and 
specifications but excluding the deviant treatments with an 
infrared laser thermometer. A total of N = 228 T „ç values 
were included in this model. 

Model 4: A more complex variant of Model 2, also con- 
sidering body size and body mass of lizards and the trial 
order of experiments as predictors. A total of N = 346 T, 
values were included in this model. 

We furthermore performed non-parametric Kruskal- 
Wallis Analyses of Variance (K-W ANOVA) with post-hoc 
tests in Statistica version 7 (Statsoft), to detect differences 
in inferred T ,¢ mong species and treatments. 


f 


Results 


After excluding experiments with outliers due to errone- 
ous measurements, data for a total of 346 experiments with 
65 individual lizards were available. For each experiment 
and individual lizard, the medians of all measurements of 
body temperature were defined as inferred I. ¿and used as 
data points for further analysis. For the photothermal- 2h 
and thermal-2h treatments which were carried out either 
in the morning or in the afternoon (see Methods), we as- 
sessed with a GLMM (Model 1) that time period had no 
significant effect on T „ (Table 1), justifying the merging 
of morning and afternoon trials in all subsequent statistics. 

To obtain a first overview of thermal preferences among 
individual lizards, and to exclude that such preferences 
would substantially influence our results, we carried out 
rank correlation analyses among the different methods 
and field body temperatures for each individual lizard, sep- 
arately by species. These analyses revealed only sporadic 
significant correlations which lost significance after Bon- 
ferroni correction (Supplementary Tables S1-83), suggest- 
ing that individual thermal preferences across treatments 
were poorly expressed. 

A GLMM (Model 2) for all 346 available Tret measure- 
ments revealed a highly significant influence of method, 
and an only weakly significant influence of species (Ta- 
ble 1). A further GLMM (Model 3) excluding the treatments 
with an infrared laser thermometer confirmed a highly sig- 
nificant effect of method whereas the effect of species was 
not significant (Table 1). A more complex GLMM (Model 
4) also considering body size and body mass of lizards and 
the trial order of experiments yielded similar results, with 
only the effects of experimental method and species being 
statistically significant (Table 1). 

Given the relatively weak or absent effect of species in 
these models, we pooled the data of all three species in 
several of the analyses reported in the following. This was 
done in an effort to increase sample size in these simpler 
statistical analyses that were carried out to better under- 
stand the differences between treatments. 

Averaged over all three species, the 2h and 10h treat- 
ments with contact thermometers (cloacal measurements 
or thermocouples attached to venter) revealed an obvious 
difference between incandescent and IR bulbs, the latter 
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Table 1. Fixed effects of various predictor variables in general 
linear mixed models (GLMMs) on Ter See Methods for detailed 
descriptions of Model 1-4. In brief, in all models T „was speci- 
fied as response variable and lizard individual as random factor. 
Model 1 was performed only for 2h treatments, and Model 3 ex- 
cluded treatments with laser thermometer measurements. Model 
4 is a variant of Model 2 with additional predictor variables. Ab- 
breviations: DF, degrees of freedom (i.e., number of fixed factors 
-1); DFDen, denominator degrees of freedom for the effect test 
(the degrees of freedom for error); F ratio, computed F ratio for 
testing that the effect is zero; Prob > F, P-value of the effect test. 


Source DF DFDen F Ratio Prob >F 
Model 1 

time (morning/afternoon) 1 58.9 0.10 0.7523 
treatment 1 48.6 114.64 <.0001 
species 2 47.7 2.09 0.1344 
Model 2 

treatment 3 278.4 62.47 <0.0001 
species 2 52.0 4.99 0.0104 
Model 3 

treatment 3 172.5 79.06 <0.0001 
species 2 56.7 1.24 0.2977 
Model 4 

treatment 5 277.4 62.24 <0.0001 
species 2 52.5 3.98 0.0247 
trial order 1 286.9 2.54 0.1123 
body mass 1 49.6 1.10 0.3002 
SVL 1 48.6 2.87 0.0967 


yielding about 5°C lower T „ç estimates (IR vs. incandes- 
cent bulb medians of 27.9°C vs. 33.9°C for the 2 h treat- 
ments; and 26.9°C vs. 32.2°C for the 10h treatments). For 
the photothermal approach, the T ver estimates obtained in 
the 2h treatment were on average 22°C higher than in the 
10h treatment (33.8°C vs. 31.6°C). The differences between 
all treatments except thermal-2h vs. thermal-1oh were 
highly significant (K-W ANOVA; P < 0.001) with pairwise 
post-hoc tests (see Fig. 1 for test details). When comparing 
results separately for each species the pattern was highly 
consistent, with IR bulb treatments yielding substantially 
lower Ta , estimates (Fig. 2). 

Measurements of field bo dy temperatures (T,) of the three 
species yielded similar median values for the two approach- 
es used (cloacal temperature, and surface temperature by in- 
frared laser thermometer): A. andreanskyi: T, avaa 28.6°C; 
Ty dorsal 28:3 ©3 P. vaucheri: Ty cloaca 271°C; Ty, dorsal 27-0 C3 S. per- 
spicillata: Ty, ayaa 272°C; T, aona, 28-3°C. We tested for possible 
differences in T, Jaa between species and found them to be 
non-significant (K-W ANOVA, P = 0.602). 

Comparing field T, asa, with T „ç inferred using the 
four contact thermometer-based methods (pooled for the 
three species; K-W ANOVA, H, 292 7 141.7121, p < 0. 001) re- 
vealed significant differences between T, and T „inferred 
with an incandescent bulb (post-hoc tests, p= i 0001 for 
both photothermal-2h and photothermal-1oh) but not be- 


tween T, and T „ç inferred using an infrared bulb (P = 1.0 
for thermal-2h, ‘and P =o. 37 for thermal-1oh) (see Fig. 2 for 
comparisons for each species). 

When scoring lizard cloacal temperatures during the 
photothermal-1oh and thermal-i1oh treatments, we includ- 
ed values from lizards that were temporarily not active, to 
obtain a daily average T „ç including also resting phases. 
To test the distribution of these phases without obvious 
thermoregulatory activity over the day, we analyzed hour- 
ly temperatures over the course of the experiment, pooled 
for all three species (Fig. 3). In the photothermal-1oh treat- 
ment we found a clear trend of decreasing lizard body tem- 
peratures over time, whereas in the thermal-1oh treatment, 
lizard body temperatures increased over time. 
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Figure 1. Measured preferred temperatures (T „p pooled for 66 
individuals of three species of lacertids (Atlantolacerta andreans- 
kyi, Podarcis vaucheri, Scelarcis perspicillata) at Oukaimeden, 
Morocco, for four different experimental procedures (lizard 
body temperature measured with contact thermometers, ei- 
ther as thermocouple attached to the venter, every minute for 
2 hours; or by cloacal measurements every hour for 10 hours; 
median body temperature of each lizard in each experiment 
used as data points for analysis). The boxplots indicate median, 
quartiles, minimum-maximum values (not considering outliers), 
and outliers as dots. Yellow colour indicates experiments where 
thermal gradients were established with incandescent bulbs, red 
colour indicates an infrared bulb as heat source and uniform cold 
lighting. Differences between methods are statistically significant 
(K-W ANOVA, H, 34 = 190.0, P < 0.0001). Differences in all pair- 
wise comparisons are statistically significant in post-hoc tests at 
P < 0.0001, except for the comparison between photothermal-2h 
and photothermal-10h (significant at P < 0.002), and thermal-2h 
and thermal-10h (not significant). 
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Treatments where lizard body temperatures in the gra- 
dient were measured using an infrared laser thermome- 
ter showed a pattern deviant from the overall one (Fig. 4). 
Here, the higher inferred A a values corresponded to 
the thermal-1ohL treatment using an IR bulb (32.6°C vs. 
25.8°C). 


Discussion 


Our study focused mainly on comparing experimental ap- 
proaches that have often been used in lacertid and other 
lizards, and to describe the amount of difference found 
among the T „estimates obtained with these approaches. 
Because the selected treatments differ at various confound- 
ed levels, i.e. heat source, method of temperature measure- 
ment, total experimental time, and interval between mea- 
surements, we cannot thoroughly discern the physical, 
behavioural and physiological causes underlying the ob- 
served differences, and in the following will develop a se- 
ries of necessarily speculative but testable hypotheses on 
these causes. However, as a main strength of our experi- 
mental approach, by subjecting a large number of individu- 
al lizards each to the various treatments in randomized tri- 
als, and by recovering similar among-treatment differences 
for three species, we provide conclusive evidence for actual 
differences between often-used approaches of T „ç estima- 
tion in lacertid and other lizards. We found (i) weak dif- 


ferences between short-term vs. long-term measurements 
under an incandescent bulb, and important and consistent 
differences between heat sources, i.e., (ii) between treat- 
ments with incandescent vs. IR bulbs. Finally, we also ob- 
served strongly different patterns when using (iii) infrared 
laser thermometers vs. contact thermometers. 

The weak differences between the photothermal-2h vs 
photothermal-ioh inferred T „values are likely explained 
by inherent differences between the treatments, and by be- 
havioural factors. We hypothesize that the 2h treatments 
targeted lizards in the initial warming-up phase and dur- 
ing the first phases of activity in the morning and after- 
noon when thermoregulation is usually a priority, given 
the known endogenous circadian activity rhythms of liz- 
ards (FOA & BERTOLUCCI 2003); our approach in these 
treatments was to even exclude lizards staying inactive for 
a long time from further analysis, similar to previous stud- 
ies (e.g., KIRCHHOF et al. 2017, GARcIA-PorTa et al. 2019). 
In contrast, the 10h treatment attempts to capture the pre- 
ferred temperature of a lizard over its entire daily cycle, 
which typically includes phases of inactivity, and phases 
when thermoregulation is not the main priority for the liz- 
ards, resulting in lower body temperatures. The 10h treat- 
ments may also exacerbate the hydric costs that can con- 
strain thermoregulation (SANNOLO & CARRETERO 2019), 
thereby favoring periods when thermoregulation is con- 
strained by other factors. This mechanism might indeed af- 
fect some of our target species, in particular Atlantolacerta 
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Figure 2. Field body temperatures (T,) and measured preferred temperatures (T „;), separately for three species of lacertids (Atlanto- 
lacerta andreanskyi, Podarcis vaucheri, Scelarcis perspicillata) at Oukaimeden, Morocco, for four different experimental procedures 
(lizard body temperature measured with contact thermometers, either as thermocouple attached to the venter, every minute for 
2 hours; or by cloacal measurements every hour for 10 hours). The boxplots indicate median, quartiles, minimum-maximum values 
(not considering outliers), and outliers as dots. Yellow colour indicates experiments where thermal gradients were established with 
incandescent bulbs, red colour indicates an infrared bulb as heat source and uniform cold lighting. 


268 


Inference of lizard preferred temperatures differ among methods 


andreanskyi which requires substantial humidity, as sug- 
gested by a previous study (S’kutFa et al. 2020). 

It is also obvious that both the 2h and the 10h treatments 
involve artificial disturbances that may influence the an- 
imal’s behaviour in different ways. On one hand, the 2h 
treatment involves attaching a thermocouple to the ani- 
mal. According to our observations, animals usually bask 
very calmly in one particular area of the gradient despite 
the attached thermocouple, but it cannot be excluded that 
in some cases this setup leads to increased activity to es- 
cape from the thermocouple annoyance. This in turn could 
lead to a higher variation of measured temperatures, and to 
higher average temperature values as the animal would not 
undergo resting activity at lower temperatures. On the oth- 
er hand, the 10h treatments involve repeated capturing and 
handling the animals (once per hour) which may lead to a 
less active behaviour over the entire phase, thus to a higher 
proportion of time spent in cooler parts of the gradient for 
resting, and consequently to lower average body temper- 
ature measurements. In agreement with the expectations 
from these differences among the approaches, the photo- 
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Figure 3. Body temperature measurements of lizards (pooled for 
all three species, Atlantolacerta andreanskyi, Podarcis vaucheri, 
Scelarcis perspicillata) measured per hour over the 10 hours of 
the (A) thermal-10h and (B) photothermal-10h experiments. 
Both correlations are statistically significant: (A) r = -0.1304, P = 
0.0027; (B) r = 0.2878; P < 0.0001. 


thermal-1oh treatment resulted in slightly but significantly 
lower T „ç estimates than the photothermal-2h treatment. 
In contrast, the magnitude of differences between heat 
sources (incandescent vs. IR) observed herein is striking 
and may in part reflect differences between thermal and 
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Figure 4. Measured preferred temperatures (T,,.,), pooled for 
three species of lacertids (Atlantolacerta andreanskyi, Podarcis 
vaucheri, Scelarcis perspicillata) at Oukaimeden, Morocco, for 
two different experimental procedures (lizard body temperature 
measured with an infrared laser thermometer every hour for 
10 hours). The boxplots indicate median, quartiles, minimum- 
maximum values (not considering outliers), and outliers as dots. 
Yellow colour indicates experiments where thermal gradients 
were established with incandescent bulbs, red colour indicates an 
infrared bulb as heat source and uniform cold lighting. The differ- 
ences are statistically significant (U-test; P < 0.001). Note that we 
here discuss this result as a possible artefact of laser thermometer 
results under the specific conditions of our indoor treatments, 
although elucidating the causes of this surprising result will re- 
quire future testing. 
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photothermal preferences. Not considering the results 
obtained with laser thermometers, T „çp values obtained 
in photothermal gradients were consistently higher than 
those in thermal gradients. We cannot readily explain 
these differences, but we speculate that this may reflect that 
heliothermic lizards such as lacertids prefer thermoregu- 
lating in bright sunlight, even if this leads to body tempera- 
tures above their purely thermal preference. On the other 
hand, in an IR-thermal gradient with uniform, constant 
artificial light, lizards may become less active than they 
would normally be, due to the lacking stimulus of a bright 
light source emitting at the same time thermal radiation. 
Therefore, in such an IR-thermal gradient, they would on 
average attain temperatures below their physiological pref- 
erence. 

Some hints on possible explanations of the observed 
pattern are provided by the trends of hourly body tempera- 
tures measured in the gradient during the 10h experimen- 
tal treatments (Fig. 3). Under such a photothermal treat- 
ment over an entire daily period, according to these obser- 
vations (data herein and numerous other studies), lizards 
in the first hours attain comparatively high body tempera- 
tures while basking, and then increasingly become active 
across the gradient, with lower body temperatures, or even 
become inactive with even lower body temperatures. These 
results may be further exacerbated by an endogenous circa- 
dian activity rhythms of lizards (FOÀ & BERTOLUCCI 2003) 
which could trigger initially higher activities in the morn- 
ing. On the contrary, under a thermal treatment, lizards in 
the beginning according to our observations (Fig. 3) oper- 
ate under relatively low temperatures and later attain high- 
er temperatures, which may be due to our use of artificial 
cold light rather than natural light with natural photope- 
riod (TosINI & AVERY 1994, 1996, TOSINI et al. 2001). How- 
ever, using purely natural light sources in our experimental 
setup was not feasible (e.g., through a large window) due to 
constraints in the available rooms at the study site. Perhaps 
also, in such a thermal treatment (especially in a setup with 
artificial light source) lizards require some time to learn 
how to efficiently up-thermoregulate in a purely thermal 
gradient, without the visual cues provided by an incandes- 
cent bulb, or by the sun in their natural environment where 
also the spatial distribution of temperatures is different. 

Further hypotheses for the factors influencing Tç as- 
sessments under different experiments of different dura- 
tion in photothermal setups can be derived from growth 
experiments in which hours of access to photothermal heat 
per day were manipulated to 6, 8, 10, and 12 hours, in Scelo- 
porus lizards found across altitudinal and latitudinal gra- 
dients (SINERVO & ADOLPH 1989, 1994). The results sug- 
gested adaptive differences in thermal “niches”, where the 
cool-adapted species may not have the capacity to use long 
thermal periods. Additional studies on S. occidentalis from 
a warm desert site and a cool site revealed a similar differ- 
ence (SINERVO 1990). Lizards from the cooler site did not 
use the full 6-hour treatment and were in refugia in 9 and 
12 hour treatments, while the populations from the warm- 
er site used the entire time in each treatment for thermo- 
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regulation. These findings, as well as those of PARANJPE et 
al. (2013) are germane in the use of short (2h) versus long 
(10h) experimental duration at revealing evolution in ther- 
moregulatory behaviour. They explain why in the studied 
lizards, short-term methods for Ts will give high values, 
while long exposure to a photothermal heat will reveal ad- 
ditional adaptations related to activity period that might 
impose severe constraints on cool-adapted species under 
climate change (SINERVO et al. 2018). We hypothesize that 
similar mechanisms also explain, at least partly, the differ- 
ences between short- and long-duration photothermal ex- 
periments in our study. 

Presumably, under field conditions, active lizards may 
largely operate under their preferred temperatures, and 
T, may, therefore, be used as a very approximate yard- 
stick to evaluate the accuracy of experimental T mea- 
surements. In our study, the temperatures from the ther- 
mal treatment (i.e., with IR bulb) were similar to field T, of 
the lizards. However, it must be considered that active liz- 
ards may have been observed during initial phases, or in- 
between episodes of basking; and indeed, in these phases 
they might be more easily spotted and captured, and will 
have lower body temperatures. Our study did not include a 
rigorous assessment of operative body temperatures in the 
field based on fully randomized sampling across the entire 
range of microhabitats used. Still, our results allow for the 
preliminary hypothesis that the thermal (IR bulb) T._ ,, esti- 
mates may more closely reflect the average body tempera- 
tures of active lizards in the field, which alternate between 
basking and being active in thermally suboptimal micro- 
habitats, e.g. when foraging. 

Lastly, the differences between T „estimates from mea- 
surements obtained by an infrared thermometer strongly 
deviated from the pattern observed with contact thermom- 
eters. An obvious and intuitive hypothesis would predict 
overall higher temperature measurements obtained from 
laser thermometers, as they measure temperature on the 
dorsal side of the lizards (exposed to the heating source), 
compared to thermocouple measurements on their belly or 
cloaca, thus on their ventral side. This prediction was con- 
firmed since under an IR bulb the laser thermometer mea- 
surements (thermaliohL) were consistently higher than 
the cloacal measurements (thermalioh). However, under 
an incandescent bulb this trend was inverted and the la- 
ser thermometer measurements (photothermaliohL) were 
consistently lower than the cloacal measurements (photo- 
thermalioh). On the contrary, when using the two types of 
thermometer with natural sunlight in the field, very simi- 
lar T, measurements were obtained. This suggests that un- 
der the experimental conditions of the thermal gradient, 
some laser thermometers may suffer from different arte- 
facts. Since we did not perform dedicated experiments on 
the underlying causes, we can here only speculate on the 
nature of these artefacts. One of these may be related to 
the substrate materials that we used to build the gradient: 
in order to have a lightweight, portable and heat-resistance 
field setup, we used silver-aluminium coated foamed plas- 
tic with a roughened wallpaper as substrate - and the re- 
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maining exposed aluminium coating at the gradient walls 
may have caused reflectance influencing the laser ther- 
mometer. Other possible biases may have been caused by 
a too wide measuring angle (HARE et al. 2007; although 
we took great care to keep this uniform across the experi- 
ments), possibly excessive measuring distance, and arte- 
facts we observed from old, partly discharged to fresh bat- 
teries. Furthermore, given that our model of infrared la- 
ser thermometer is recommended only for altitudes up to 
2000 m above sea level (and our study site was located at 
2600 m above sea level), we cannot exclude that using dif- 
ferent models of laser thermometer these factors and arte- 
facts can be attenuated or avoided. With the data at hand, 
it seems impossible to clearly understand what has caused 
the counter-intuitive laser thermometer measurements of 
lizard body temperature. Given these erratic results, our 
study calls for caution when using laser thermometers un- 
der experimental conditions rather than natural sunlight. 
Overall, our comparisons highlight that experimental 
approaches to determine T,,,, can yield substantially differ- 
ent results. First and foremost, our experiments confirmed 
the existence of important differences between laser and 
contact thermometers (CARRETERO 2012), and call for ex- 
treme caution and need for careful exclusion of possible 
measuring artefacts when using the former instruments 
especially in non-natural setups. Second, we obtained sur- 
prisingly high differences in median values between setups 
with incandescent vs. IR bulbs (photothermal vs. thermal), 
with average differences of about 5°C in inferred T „p This 
differs from empirical differences summarized by GARCIA- 
Porta et al. (2019) who found an average difference be- 
tween these methods of only 1°C, and applied this value as 
correction factor in order to include in their meta-analysis 
a few crucial species for which only data from a 10h-ther- 
mal setup were available. As discussed above, this may be 
due to the use of non-natural light in our thermal setups, 
and highlights the need to consider lighting as important 
factor when interpreting the results from thermal gradient 
Tę assessments. Thirdly, it is encouraging that T, data 
inferred from different photothermal treatments (2h vs. 
10h) under incandescent bulbs, and obtained with contact 
thermometers, resulted to be roughly comparable, at least 
in the three sympatric lacertids studied herein, with a sug- 
gested correction of +2.2°C to make the photothermal-1oh 
comparable to the photothermal-2h treatment. Yet, given 
the known inter- and intraspecific differences in thermo- 
regulatory period (see above; SINERVO & ADOLPH 1989, 
1994, SINERVO 1990) it is unlikely that such a correction 
factor can ever be generally valid for all lacertid species, let 
alone for other lizards. As previously recommended for la- 
ser and contact thermometers (CARRETERO 2012) and now 
extended to the other methods, where detailed and thor- 
ough comparisons are needed, a rigorous consideration of 
the optimal experimental setup and its consistent use will 
remain necessary. Lastly, we envisage that future field and 
experimental studies on lizard thermoregulation will in- 
creasingly make use of thermographic cameras which al- 
low visualizing temperatures even in different lizard body 


parts, without disturbance, but require careful calibration 
of the equipment for the target species (BARROSO et al. 
2016). 
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Supplementary data 


The following data are available online: 


Supplementary Table S1. Matrix of Spearman rank correlation 
coefficients for median T „values obtained from different meas- 
uring methods, and field t body temperatures (T,), for 30 male 
lizard individuals of Atlantolacerta andreanskyi at Oukaimeden. 


Supplementary Table $2. Matrix of Spearman rank correlation 
coefficients for median T „values obtained from different meas- 
uring methods, and field | body temperatures (T,), for 23 male 
lizard individuals of Podarcis vaucheri at Oukaimeden. 


Supplementary Table S3. Matrix of Spearman rank correlation 
coefficients for median T „values obtained from different meas- 
uring methods, and field | body temperatures (T,), for 12 male 
lizard individuals of Podarcis vaucheri at Oukaimeden. 
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